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The low energy behavior of the 5 = 1/2 antiferromagnetic Heisenberg chains with precious 
mean quasiperiodic exchange modulation is studied by the density matrix renormalization group 
method. Based on the scaling behavior of the energy gap distribution, it is found that the ground 
state of this model belongs to the universality class different from that of the XY chain for which 
the precious mean exchange modulation is marginal. This result is consistent with the recent 
bosonization analysis of Vidal et al.™ 
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In the recent studies of quantum many body problem, 
the low energy properties of the quantum spin systems 
with modulated spatial structure have been attracting 
a wide interest. Although the periodic chains and the 
random chains are studied in detail, the quasiperiodic 
chains, which has the intermediate character between 
the regular and random chains, are less studied except 
for the XY-casa whiiitL is equivalent to the spinless free 
fermion chainsoBuurcP In the fermionic language, the 
Ising component of the exchange coupling corresponds 
to the fermion-fermion interaction leading to the strong 
correlation effect which is the most important subject of 
the recent condensed matter physics. 

Although the XY chain can be mapped onto the free 
fermion chain, the problem is not trivial on the quasiperi- 
odic lattice- rjjEpr the Fibonacci lattice, Kohmoto and 
coworkersoBtlcP clarified the Cantor set structure of the 
single particle spectrum and the wave function by means 
of the renormalization group (RG) method. Especially 
the dynamical exponents arc found analytically at the 
band center and the band edge. Recently, this approach 
has been extended to include other types of quasiperiodic 
lattices and the anisotropy between x and y component 
of exchange couplings© It should be noted that the criti- 
cality of the Fibonacci XY chain stems from the marginal 
nature of the Fibonacci and other precious mean aperi- 
odicity in this model. For the relevant aperiodicity, more 
singular behavior with divergent dynamical exponent is 
realized even for the XY chain Jy 

Although these works discovered the beautiful mathe- 
matical structure of quasiperiodic chains, almost no at- 
tempt to include the interactioiLeffect is done so far ex- 
cept for the mean field approachauP and recent bosoniza- 
tion approach!! In the present work, we employ the der 
sity matrix renormalization group (DMRG) methocoL 
to take full account of the correlation effect in the S 



1/2 precious mean antiferromagnetic Heisenberg chains 
which include the Fibonacci chains. 
Our Hamiltionan is given by, 

N—l 

H=J2 2J <*iSiSi+i, (J at >0), (1) 
»=i 

where Si's are the spin 1/2 operators and the open 
boundary condition is assumed. The exchange couplings 
•Toi's (= Jk or Jg) follow the precious mean sequence 
generated by the substitution rule, 



A — » A B. B 



A. 



(2) 



The cases k = 1 and k — 2 correspond to the Fibonacci 
(golden mean) and silver mean chains, respectively. In 
the following, we take J a = 1 to fix the energy unit. For 
finite N, we consider all possible (N — l)-membered sub- 
sequences of the infinite precious mean chain and investi- 
gate the energy gap distribution among them. It should 
be noted that the number of the (N — l)-membered sub- 
sequence is equal to N£3 

In the XY case, the precious mean aperiodicity in the 
exchange coupling is marginal and the energy gap A 
scales with the system size as A ~ N~ z where z is the 
dynamical exponent o Therefore the gap distribution 
function scales as 



P{A)dA = N z f(AN z )dA. 



(3) 



Consequently, the average and flucutation of ln(l/A) 
scales as, 



< ln(l/A) > ~Ci - zlnN, 



ffpn(A)] = V< (MA)- < m(A) >f >, 
C2 — C\ = const., 



(4) 
(5) 
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where 

/oc 
t n g(t)dt, g(t) = fie"). 
-oc 

On the other hand, for the XY chain with relevant ea 
change aperiodicity, the gap scales as ln(l/A) ~ N^IL 
and the gap distribution function scales as, 



P(lnA)dlnA 
which gives 



7V~"/(iV- w lnA)dlnA, 



< ln(l/A) > ~ DiN u 



where 



a[ln(A)] 



D n = 



D 2 - D\N U 



x n f(x)dx. 



(6) 



(7) 
(8) 



It should be remarked that a tends to a constant value for 
the marginal aperiodicity while it grows with the same 
exponent as < ln(l/A) > for the relevant aperiodicity. 
This type of behavior with co = 1/2 is observed also in 
the random singlet phase of the_|S-.= 1/2 random anti- 
ferromagnetic Heisenberg chain.BliJ 

The interacting spinless fermion chain with Fibonacci 
potential has been studied by Vidal et al£P by means of 
the bosonization technique. This model can be mapped 
onto the XXZ chain in the Fibonacci magnetic field by 
the Jordan- Wigner transformation. In the following, we 
use the spin chain terminology. After bosonization, this 
model can be described by the boson Hamiltonian 



Vidal et alxS derived the RG equation within the weak 
coupling approximation. Based on the numerical solu- 
tion of the RG equation, they have shown that the Fi- 
bonacci modulation term becomes irrelevant or relevant 
according as K > K c or K < K c where K c ~ 2 in the 
absence of the uniform magnetic field. For K > K c , 
the ground state is the usual Luttinger liquid, while for 
K < K c , the ground state is renormalized to the strong 
coupling regime and the weak coupling theory cannnot 
predict the ground states properties. For K — K c , which 
corresponds to the XY model, the Fibonacci modulation 
becomes marginal and the ground state is critical. This 
is consistent with the well-known case of the free spinless 
fermions in the Fibonacci potential. On the other hand, 
the isotropic case (K = 1) is renormalized to the strong 
coupling regime. 




Ho — — I dx 



Tjh 



2tt 



2na 



H = H + H^, (9) 



{uK ){^f + [-){d^) 2 , (10) 



dx W(x) cos \2k F x + V2(j)(x)] , (11) 



where <j> is the boson field, II, the field conjugate to tj>, 
a, the ultraviolet cut-off, u, the spin wave velocity, hp, 
the fcrmi wave number of the spinless fermions and K, 
the Luttinger liquid parameter. The function W(x) rep- 
resents the spatially varying magnetic field. The case 
K = 1 corresponds to the SU(2) invariant isotropic 
Heisenberg chain and K — 2to the XY chain. (Note that 
our definition of K differs from that of ref. ^ by a factor 
of 2.) For the Fibonacci type modulation, the function 
W{x) is defined via its Fourier components given in ref. 



Using the standard bosonization scheme, the spatial 
modulation of the exchange coupling is similarly ex- 
pressed as, 



Hw = 7T~ I dxW(x) sin 
2ira 



2k F x + V2(f>(x) , (12) 



which coincides with eq. ([11]) by the shift of the origin of 
(j>. Therefore the conclusion obtained by Vidal et al. also 
holds for the case of Fibonacci type exchange modulation 
given by the Hamiltonian (Q). 



Fig. 1. The jV-dependence of < ln(l/A) > and <r[lnA] for the 
Fibonacci XY chain by the exact diagonalization method. The 
length of the horizontal bar is 31n((l + v / 5)/2). In this and 
following figures 3, 5 and 6, r = Max{ Jg/J^, ^a/^b} an< ^ * ne 
filled (open) symbols represent the case Jg > J\{J-q < J^)- 
The filled symbols for a almost overlap with open symbols. 




Fig. 2. The dynamical exponent z for the XY chain. The filled 
circle represent the present numerical caloHWtion and solid line 
is the analytical results by Kohmoto et aZ.ErU 



In the following, we present the numerical results on 
the energy gap distribution. In the XY chain, the energy 
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spectrum can be calculated numerically by the diagonal- 
ization of N x N matrices. The average < ln(l/A) > and 
fluctuation <j[lnA] of the logarithm of the energy gap for 
all possible (N — l)-membered subsequence is calculated 
for 14 < N < 378 and various values of Jg/ 'Ja between 
1/5 and 5. These are plotted against h\N in Fig. 1. The 
linearity of < ln(l/A) > to hiN is fairly good and the 
flucutation er also tends to a constant value except for the 
oscillation with period 31n{(l + \ZKl/2} which is inherent 
to the Fibonacci chain spectrumJlStP The facter 3 comes 
from the fact that the single step of the RG transforma- 
tion for the precious mean chains with odd k corresponds 
to three inflation steps£P We further calculated the dy- 
namical exponent z by fitting ln(l/A)-ln./V curve by a 
straight line. The obtained values of z are shown in Fig. 
2 by open circles. The solid line is the^xact expression 
obtained by Kohmoto and coworkersotP Again we find 
a good agreement. This result also confirms that the en- 
ergy gap distribution of the finite length subsequences of 
the Fibonacci chain correctly reflects scaling properties 
of the ground state in the thermodynamic limit. 




Fig. 3. The ^-dependence of < ln(l/A) > and <x[lnA] for the Fi- 
bonacci Heisenberg chain by the DMRG method plotted against 
InTV. 



The energy gap distribution for the Fibonacci Heisen- 
berg chain is calculated by the DMRG method for 14 < 
N < 234 and various values of Jg/JA between 2/3 
and 3/2 using the algorithm developed for the random 
chains .EP The number m of the states kept at each step of 
DMRG iteration is 60 and only infinite size iterations are 
carried out. We have checked that m-dependence is neg- 
ligible by increasing m up to 80 for Jb/Ja = 3/2 which 
is the most dangerous case studied here. For the further 
check of the accuracy of the DMRG scheme, we have also 
calculated the energy spectrum of the XY chain with 
14 < N < 234 using DMRG and found that the results 
coincide with the exact diagonalization results within the 
size of the symbols of Fig. 1. 

Figure 3 shows the average < ln(l/A) > and fluctua- 
tion cr[A] plotted against IniV for the Fibonacci Heisen- 
berg chain. The curves of < ln(l/A) > show an evident 



upturn as N increases. Further, the fluctuation a does 
not tend to a constant value. We have fitted the data 
for < ln(l/A) > by the power law N u . Such behav- 
ior is, expected for the relevant aperiodicity in the XY 
case.cP The exponent uj turned out to be non-universal 
depending on the ratio Jg/J/^ as to = lu(J-q/ JpJ. We 
have further assumed that w(Jg/ Ja) = lo(Jj^/ Jg) be- 
cause the RG equation of ref. [7] is invariant under the 
exchange Ja <-> Jg. It should be also noted that the 
dynamical exponent for the precious mean XY chain is 
also invariant under the exchange Ja <-> Jg . In addition, 
the numerically obtained values of u> for the silver mean 
chains also satisfiy the relation w(Jg/ Jj^) ~ w(J^/Jg) 
as explained below. So we assume that this relation also 
holds for the Fibonacci Heisenberg chains and use the 
average of the numerically obtained values of w(Jg/JA ) 
and cj( J a / Jg) as to. The values of a; are depicted against 
lnr (r = Max(Jg/J^, J^/Jg)) in Fig. 4 by filled sym- 
bols. The error bars in Fig. 4 are estimated from the 
difference between to(J-o/J\) and to{J\/J-n). 
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Fig. 4. The exponent u for the Fibonacci and silver mean Heisen- 
berg chains. 



Using thus obtained values of to, < ln(l/A) > and 
a - [In A] are plotted against TV" in Fig. 5. It is clearly 
seen that both < ln(l/A) > and cr[lnA] grow linearly 
with . The same analysis is made for the silver mean 
chain in Fig. 6 for 18 < N < 240. In this case, the oscil- 
lation period in InN is \n(l + V2), because the single step 
of the RG transformation for the precious mean chains 
with even k corresponds to a single inflation stepS* In 
this case, the numerically obtained values of to satisfy 
w(J^/Jg) ~ (^(J-q/Jj^) with better accuracy than the 
Fibonacci case. The values of to are plotted against lnr 
in Fig. 4 by open symbols. 

The above behavior of the energy gap clearly suggests 
that the ground state of the precious mean Heisenberg 
antiferromagnetic chain belong to the universality class 
different from that of the XY chain for which the precious 
mean aperiodicity is marginal. This result is consistent 
with the bosonization analysis of Vidal et al.EP Although 
these authors predict nothing about the ground state 
in the strong coupling regime, our numerical calculation 
shows that the gap distribution is characterized by the 
scaling form (^|) which is similar to that of the XY chain 
with relevant aperiodic modulation.EP It should be re- 
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Fig. 5. The ^-dependence of < ln(l/A) > and <r[lnA] for the Fi- 
bonacci Heisenberg chain by the DMRG method plotted against 




Fig. 6. The JV-dependencc of < ln(l/A) > and <r[lnA] for the 
silver mean Heisenberg chain by the DMRG method plotted 
against . 



marked in the random exhange case, the S — 1/2 XXZ 
model (1 < K < 2) belong to the stroag^Ottpling regime 
because the critical value of K is 3J3't3'EJ> For these 
models, the ground state is known to be the random 
singlet state which is also characterizcdJajy the gap dis- 
tribution of the type (||) with u> = 1/2H Presumably, 
the gap distribution of type (^|) is the generic nature of 
the aperiodic exchange spin chains in the strong coupling 
regime. 

To summarize, based on the DMRG calculation of the 
gap distribution, it is found that the logarithm of the 
energy gap of the precious mean Heisenberg chains scales 
with a nonuniversal power of the system size as eq. ( |q) . 
This is distinct from the critical gap distribution eq. (H) 
in the XY case. Taking_±h£. results for other kinds of 
aperiodic spin chainsQBlHrEj) into account, we expect 
that this is the characteristic property of the aperiodic 
spin chains in the strong coupling regime. 



In this paper, we concentrated on the energy gap dis- 
tribution in the absence of the magnetic field. The ef- 
fect of the magnetic field is interesting from two different 
points of view. First, in the presence of the uniform mag- 
netic field, the multifractal Cantor-set structure of the 
single particle, .excitation spectrum of the free spinlcss 
fermion chairlltLj) manifests itself as the devil's staircase 
structure of the magnetization curve of the XY chain in 
the spin language. This_caii be regarded as the magneti- 
zation plateau problcmOEZP in the quasiperiodic spatial 
structure. It is worth investigating if such structure sur- 
vives in the Heisenberg case which has stronger quantum 
fluctuation than the XY case. 

Another interesting problem is the effect of the pre- 
cious mean modulation of magnetic field. For the ran- 
dom field XY chain, the band center state is known to 
be localized while it is the random singlet state with di- 
vergent spin correlation length for the random exchange 
XY chain, although the randomness is relevant in both 



This difference comes from the perfect s 



cases 

inversion symmetry of the random exchange problem. I 
From this point of view, effect of the precious mean mod- 
ulation of the magnetic field might be different from that 
of the exchange coupling in the Heisenberg chain. This 
problem is left for future studies. 
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